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I. INTRODUCTION 

In recent years, a significant research effort has been dedicated to the development of new flow energy-harvesting 
systems and the characterization of their efficiency. Many of such systems are based on flow-induced vibrations of 
structures, such as vortex-induced vibrations m, flutter of wing profiles HH, of cylinders in axial flows BIS], 
or flags EHia. The present work focuses on the latter, and considers a flag covered with piezoelectric patches that 
convert mechanical deformation into electric charge transfer inside an output circuit. 

In this article the term flag refers in fact to a clamped-free plate in an axial flow. It is well-established that 
such structure can undergo strong self-sustained oscillations once a critical flow velocity is overpassed (see the recent 
review by Ref. |13] and references therein). These oscillations originate from an instability of the equilibrium position, 
resulting from the interaction of the flow forces with the solid’s inertia and rigidity. This instability can be predicted 
from the linearized dynamics equations and different methods and models have been proposed depending on the aspect 
ratio of the flag, ranging from two-dimensional models (large span, [T31[T5]) to slender body models (short span, [IBjb 
The general case of a three dimensional fluid flow around a rectangular plate showed the limits of applicability of 
both models HZIIIH]. Linear analysis is however not able to predict the amplitude of flapping in the saturated regime, 
hence it is necessary to consider the nonlinear dynamics of the flag to be able to determine how much energy can 
be harvested. Beyond full numerical simulations of the coupled fluid-solid problems, simplified fluid-solid nonlinear 
models have been proposed in the limit of two-dimensional flows [El [JQ] or slender body problems mm- 

Piezoelectric patches have traditionally been studied to couple structural mechanics to electrical circuits in the 
context of passive damping of structural vibrations |22]j of active control of vibrations [23] and of course energy 
harvesting [24ll27j . The use of a piezoelectric flag to convert the kinetic energy of a flow was initially introduced by 
Ref. [T] , to produce energy from the flapping of a flexible membrane forced by the wake of an upstream obstacle. In the 
last five years, several experimental, numerical and theoretical studies have been conducted to study the possibility 
to exploit the flutter instability itself using piezoelectric plates [El dHUHl [29] . In these works, the plate was either 
entirely or partially covered by one piezoelectric pair [281129j . or completely covered by a large amount of small 
piezoelectric elements, so that a continuous model for a homogenous plate can be considered [E1[II|. In the latter 
case and for purely resistive output circuits, it was shown that the maximum efficiency is obtained for large mass 
ratio when the timescale of the electrical circuit is equal to the fluid-solid instability timescale m- For practical 
applications, however, only a finite number of piezoelectric elements can realistically be used. The goal of the present 
article is therefore to provide some understanding on the impact of such a discrete coverage on the dynamics and 
performance of the system, a critical question that has remained so far unanswered. 

In contrast with simplified representations of the energy harvesting process (e.g. pure damping, [30|), such an 
approach provides a fully-coupled description of the fluid-solid-electric problem, taking into account both the transfer 
of energy from the solid motion to the output circuit and the feedback effect of the piezoelectric patches on the 
flag’s dynamics. The nonlinear dynamics of the flag (even in the absence of any piezoelectric coupling) exhibit strong 
inhomogeneities of the local deformation and motion along the plate: near the clamped leading edge, the flag has little 
displacement but significant curvature, while near the free trailing edge, the flag has zero curvature but undergoes 
large displacements. Local stretching and compression of the piezoelectric patches resulting from the flag deformation 
are responsible for the charge transfer in the output circuit; the location along the structure of the device converting 
mechanical to electrical energy is therefore expected to strongly affect the efficiency of the system. A similar question 
was recently addressed in Refs. [Hi in the case of a pure damping model. 

The present work focuses on the optimization of the distribution of a small set of piezoelectric elements (one, two 
or three pairs) positioned on a plate in an axial flow. In section |nj the physical model and equations are presented 
for the fluid-solid-electric system with a finite number of discrete piezoelectric patches. A method of simulation based 
on the weakly nonlinear form of these equations is presented that is appropriate to obtain the nonlinear dynamics 
near the instability threshold. The influence of the number of piezoelectric elements on the stability of the system 
is investigated in section [ml In section EYI the impact of this arrangement on the nonlinear dynamics and energy 
harvesting efficiency is studied using experiments and numerical simulations. Finally, a parametric numerical study 
is performed to determine the optimal arrangements in the case of one, two or three piezoelectric elements. 


II. PHYSICAL MODEL 

The system considered in this work consists of an elastic flexible plate of length L and width H immersed in an axial 
flow of uniform density p and velocity Uoo ■ The plate is inextensible and clamped at its leading edge. For simplicity only 
purely two-dimensional motions of the plate are considered. The plate’s geometry is therefore completely determined 
by the local orientation of the plate with respect to the flow, 9{S,T), where S is the curvilinear coordinate along the 
plate and T is time. A finite number Np of piezoelectric electrodes of length Lpi and width H are attached on each 
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FIG. 1. (a) Two dimensional flapping of a flexible plate covered by pairs of piezoelectric patches, (b) Flexible plate flapping in 
uniform axial flow. 


side of the plate. For each piezoelectric patch pair, the negative electrodes are shunted across the plate while their 
positive electrodes are connected to an external circuit. The charge (per unit length in the streamwise direction) in 
each piezo pair is given by m- 


Q^ = ^m%+cv,. ( 1 ) 

Ljpi I 

with Vi the voltage between the positive electrodes of the fth piezo pair whose left and right edges are positioned 
at S~ and respectively. C = CijLpi is the equivalent capacitance per unit length of the piezo pair and x is the 
mechanical/piezoelectrical conversion factor of the piezoelectric material. Considering that each electrode is connected 
to a purely resistive circuit, Ohm’s law leads to: 


+GV -0 


( 2 ) 


where G = Gi/Lpi = 1/RiLpi is the conductivity per unit length of the harvesting circuit, with Ri the circuit’s 
resistance. 

The voltage V) between the positive electrodes also generates an internal torque in the piezoelectric patches, and 
thus on the flexible plate, so that the total internal torque in the plate and piezoelectric patches assembly is obtained 
as 
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( 3 ) 


where B is the flexural rigidity, which depends on the respective thickness. Young’s modulus and Poisson’s ratio of 
the two constitutive materials [101. equation ([^, Fi{S) is the polarization function of the ith patch. In the present 
approach, Fi{S) = Hs{S — S~) — Hs{S — ), with Hg the Heaviside step function. The conservation of momentum 

then leads to m- 


VF,j nj - Pn, (4) 

where X(iS', T) is the local plate position and ^ is the mass per unit length of the assembly. Ft is the internal tension 
in the plate that enforces inextensibility and P is the pressure force density exerted by the surrounding fluid over the 
plate. In Eq. Q, n and r are the unit normal and tangent vectors to the plate (see Figure [^. Following previous 
works on flapping flag modeling P is decomposed into a reactive part Preac accounting for potential flow effects 

(e.g. added mass), corrected by a resistive part Pres accounting for the form drag on the plate resulting from lateral 
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flow separation: 


P = Pr. 


-fres — rn 


apH"^ - {wuy + + ^pCdH\l 


(5) 


where rUa is the added mass coefficient of the transverse section {rua = 7 r /4 for a rectangular plate |3I]), and Cd = 1.8 
the drag coefficient |32j . w and u are respectively the normal and longitudinal components of the plate’s velocity 
relative to the fluid flow, such that X — Uoo^x = ut + icn, with the unitary vector in the x direction. 

In the present approach, it is assumed that the electrode does not affect the density and bending rigidity of the 
sandwich plate that are considered homogeneous along the assembly (see [33] for a study on the passive damping of 
vibrating beams including this effect). Also, in the approach followed in this work, P^es accounts for the effect of the 
lateral flow detachment in the form of a drag. Flow detachment from the trailing edge and wake dynamics are only 
accounted through the reactive term which models the advection of added fluid momentum (see Ref. |34] for more 
details). Other models exist that consider the influence of an unsteady wake [THl 135] . 


A. Non-dimensional equations 


Using L, L/Uao, Uao\/~pJc and Uaoyfpc as characteristic length, time, voltage and charge respectively, equations 
0-0 become in non-dimensional form: 


<li = V^ + 77 ^ [ 6 ']^*-, 


U*7i 


’w+ 


. 1 92 /QQ 


dP ~ 5s U*'^ ds'^ V5s'^' Pres+M H WaPreac) n 


+ ^ H - Si ) - - 4 )] n. 


( 6 ) 

( 7 ) 

( 8 ) 


with clamped-free boundary conditions. The problem is characterized by the following non-dimensional parameters: 
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- the non dimensional velocity U* = LUao\ —i 

B 

- the electro-mechanical coupling factor a = 


X 




- the fluid-solid inertial ratio M* = 


pLH 


H 


- the non dimensional plate span H* = —, 

L 

- the non dimensional electrode length 7 ^ = —^, 

Ij 

- the tuning coefficient of the fluid-solid and electric systems f3 = 


U^C 
LG ■ 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


B. Weakly non-linear form of the eqnations 


Equation ([^ is projected onto the x and y directions in order to obtain two equations for x{s, t) and y{s, t) respectively. 
The horizontal projection is used to eliminate the tension term /,- from the y projection. Finally, x and its derivatives 
are eliminated using the inextensibility condition. Keeping terms up to 0{y^) one obtains a weakly non-linear equation 
for y{s,t) (see also Ref. [3T] for the case of a non-piezoelectric plate): 


L{y) + fmiy) + 


1 

{7*2 


/b(2/) - 


^/x(y> + M*f^es{y) + M*H*mafrea.c(,y) = 0 , 


(15) 
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where 


L{y) = y+ - sf)) + M*H* {y" + 2y' + j/), 


a 

U 


(16) 


pH ps 


{y'^ + y'y')dsds', 


fm{y) = yi + y'y')ds - y" 

Jo Js Jo 

/B(y) = wy" + y'V"^ + ^/"^ 

f^{y,vi..VNp) = \y'‘^'^Vi{5'{s - s~)- 5\s - sf)) + y'y'''^v,{6{s - s~) - 5(s - sf)), 


(17) 

(18) 
(19) 


/res(j/) = 2^D\y' + 2)1 (y' + y), 

/reac(y) = “^2/"^ ^ + v '“ 3y"y'y - 2^'/?) - ^y"y^ J ^ 

+ 2(y" + y') / y'y'ds-y" f y'{y" + 2y' + y)ds. 


( 20 ) 


( 21 ) 


Equation (16) corresponds to the linearized dynamics while Equations @-(1111) correspond to nonlinearities related 
to inertia, stiffness, piezoelectric coupling, resistive and reactive flow effects, respectively. 


C. Numerical method 


A Galerkin decomposition is used to solve Equation (151: the vertical displacement y is expanded as a superposition 
of clamped-free beam eigenmodes y(s), 


y{s,t) = ^ Ap(f)(/)p(s). 


( 22 ) 


Next, equation (151 is projected on the same set of eigenmodes. Details of the Galerkin projection can be found in 
[Aj After truncation to N linear modes, the resulting coupled system of equations is integrated numerically using a 
semi-implicit step-adaptive fourth order Runge-Kutta method. 


D. Harvested energy and harvesting efficiency 

Once the dynamics of the piezoelectric flag is obtained using the method above and a limit-cycle oscillation is reached, 
the total energy harvested by the system is computed as the energy dissipated in all the resistive circuits. In dimen¬ 
sional form, this reads as 


Pe = 



(23) 


with (.) the temporal average over one oscillation period. The harvesting efficiency y = Pe/Pf is then computed as 
the fraction of the fluid kinetic energy flux Pf through the vertical cross-section occupied by the flag that is actually 
harvested and dissipated in the output circuits. Using the non-dimensional variables defined above, this can be written 
as 


1 




Tj = 


M*(3A* 




\i^l 


(24) 


with A* = AjL and A the flapping amplitude of the flag. 
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FIG. 2. (a) : Evolution of the stability threshold as a function of M* for different values of the coupling factor a, /3 = 1 = cst, 
Np = 100. (b) : Impact of the number of electrodes on the critical stability threshold. Values of U* are normalized using the 
critical value corresponding to Np = 100. 


III. FLAPPING INSTABILITY : LINEAR ANALYSIS 

As a first validation of our model, we investigate the effect of the piezoelectric arrangement of the flag on its stability 
for a continuous coverage = sf), and a large number of piezoelectric patches. Beyond a critical value U* of the 
flow velocity the system becomes unstable to flutter [T31. Figure (left) shows the critical velocity U* as a function 
of the mass ratio M* for different values of the coupling coefficient and a large number of small piezoelectric patches 
{Np = 100). Without piezoelectric coupling, U* is a decreasing function of M*. Cusps in the curve represent a switch 
in the mode becoming unstable at the lowest flow velocity (see for example m)- These results match that obtained 
in Ref. m using a continuous model, thus indicating that for a large amount of piezoelectric patches, the present 
discrete model indeed converges to the continuous framework. 

The electro-mechanical coupling also clearly influences the stability threshold of the flag. In the coupled case 
{a 7 ^ 0) and for small values of M*, the damping induced by the electro-mechanical coupling increases the critical 
velocity U*. This effect is more important as the coupling factor a is increased. However, for high values of M*, the 
electro-mechanical coupling has a destabilizing effect, which was interpreted in Ref. m in terms of negative energy 
waves present on the non-dissipative flag [361138] . 

Several previous studies have focused on the limit of the continuous coverage of the plate by infinitesimal patches 
{Np ^ 1, [ini dH IM])- The discrete approach considered here allows us to investigate the convergence of the results 
with Np and the applicability of the continuous approach to experimental situations where a finite and often small 
number of patch pairs is considered. Figure (right) indeed shows a rapid convergence of the critical velocity for 
Np > 12, beyond which the results of the continuous model are recovered [iniiiiiEn]. Moreover, it also shows that 
using fewer patches can significantly modify the stability threshold, positively or negatively, which motivates the 
present study on a finite number of piezoelectric pairs. 

The rapid convergence of the results beyond Np « 10 can be understood as follows: for the values of the mass 
ratio considered, the dominant unstable mode is of low order resulting in typical wavelengths greater than half a flag 
length. For a large enough number of patches, the patch length is small enough that the entire mode structure is well 
captured and the discrete nature of the piezoelectric coverage has no influence. 


IV. OPTIMIZATION OF THE ELECTRODES’ POSITION 

During its self-sustained flapping motion, the flag undergoes complex and non-uniform deformations. Electric 
charge transfer between the patches’ electrodes is intimately linked to the local deformation of the flag. An important 
challenge for maximizing the amount of energy harvested by such a device therefore lies in the correct positioning of 
the electrodes. When the electromechanical coupling is large, the piezoelectric patches can significantly modify the 
flag’s dynamics resulting in a non-trivial optimization process. 
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In order to address this question, we first present a set of preliminary experiments exploring the impact of electrode 
size and position on the global energy harvesting efficiency of the system. Experimental results are then compared to 
numerical simulations based on the physical model presented in section|n] In the last part of this section, a parametric 
study is performed based on numerical simulations using more complex geometrical configurations, to identify optimal 
positioning of the electrodes. 


A. Experimental study of the role of electrode size and position 

The aim of the experiments described below is to obtain some preliminary evidence of the impact of size and position 
of the electrodes on the harvested energy. For this purpose we have focused on a simple configuration consisting of 
a flexible piezo-electric plate covered by a single pair of electrodes of variable length and position, immersed in an 
axial wind flow of speed (figure [^a). Two different experimental configurations are considered (see figure [^b for 

a graphical description of the experiments): 

— Experiment 1: the length of the electrode Lp is varied while one of its extremities is fixed at the clamped edge 
of the plate, 

— Experiment 2: the length of the electrode is also varied but with one if its extremities fixed to the free edge of 
the plate. 


(a) 




H 



(c) 

BOPP layers 



(d) 





double sided 
scotch tape 


piezoelectric 

film 


FIG. 3. (a) Description of the experimental setup, (b) Graphic representation of the electrode size and distribution for 

experiments 1 and 2. (c) Schematic of the multilayer composition of the piezoelectric flag, (d) Connection diagram of the 
harvesting circuit. 


In these experiments, the plates include a center core consisting of a double sided adhesive tape of length L = 12 
cm, width H = 2.5 cm and thickness hg = 100 /rm. Both sides of the tape are partially covered by a piezoelectric 
PVDF film of the length L and width Hp = 0.3H, with its negative polarity facing inwards. The thickness of the 
piezoelectric films is hp = 40 /rm for experiment 1 and hp = 50 /im for experiment 2. The remaining free surface at 
both sides of the scotch tape is covered by a layer of Innova biaxially oriented polypropylene (BOPP) film of exactly 
the same thickness as that of the piezoelectric film. Figure |^c shows a schematic view of the plate composition. The 



































































FIG. 4. Top left: Image sequence of the flag deformation over one oscillation period. Top right : Voltage signal generated by 
the piezoelectric patches during a typical experiment. Bottom : Normalized dissipated power Pjv as a function of the tuning 
parameter /3. Data corresponds to experiment 1, different symbols correspond to different values of Lp. The normalized power 
Pn is calculated as Pn = P/Pmax for each given value of Lp. The solid line corresponds to a numerical simulation using Lp = 1. 


piezoelectric PVDF films are covered with an external Cr/Au which serves as the electrode. In both experiments this 
Cr/Au layer is sectioned in order to adjust the desired value of <5 and Lp. The objective of covering only one third of 
the width of the plate with piezoelectric film is to reduce the impact of the sectioning process on the global rigidity 
of the plate. 

For each experiment and a given electrode length and position, the experimental procedure is as follows : the 
piezoelectric plate is placed in a wind tunnel of rectangular transversal section 10 cm wide and 5 cm high with 
transparent walls that allow visual access from the outside m- The negative electrodes of the piezoelectric films are 
shunted while a resistance R is connected between the positive electrodes (see figure |^d). Once the system is set we 
proceed to generate a wind flow inside the tunnel. Initially at rest, the wind velocity is gradually increased until the 
plate starts to flap in a self sustained oscillatory regime (upper left image in figure]^, at which point the voltage 
Vo generated along the resistance R is measured (figure upper right). The power dissipated in the resistance is 
then calculated as the temporal average Pe =< VolR>, for a 5 second time frame. The resistance R is varied in 
order to find the value of j3 that maximizes the dissipated power. Indeed, it is expected that the harvesting efficiency 
or damping presents a maximum between the extreme values /3 = 0 (short circuit condition) and j3 = oo (open 
circuit condition). The value of (3 maximizing the efficiency corresponds to a tuning of the flapping frequency to the 
characteristic timescale of the RC circuit consisting of the piezoelectric capacitance and output resistor pT] . Figure]^ 
(bottom) shows a good agreement between experiments and simulations obtained from the model described in section 
[nl The discrepancies observed between the experimental and numerical results are mainly due to small differences in 
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FIG. 5. (a) : Harvesting efficiency as a fnnction of the electrode length Lp for experiment 1. Experimental parameters are : 
U* = 14.6, M* = 0.68, H* — 0.3, /3 = 0.8 and a = 0.073. (b) : Harvesting efficiency as a function of the electrode length 
Lp for experiment 2. Experimental parameters are : U* = 15.7, M* = 0.61, H* — 0.3, P = 0.8 and a = 0.062. Error bars 
correspond to the error propagation of the standard deviation of Vo used to calculate the dissipated power Pe =< Vo IR > 
during a 5 second time frame. Solid line corresponds to numerical simulations. 


the flapping frequency of the flag, mainly caused by the error in the estimated value of the plate rigidity B used in 
the simulations. It is important to remark that the optimal /3 value remains fairly constant regardless of the size and 
position of the electrodes, a direct consequence of the flapping frequency remaining unchanged when changing the 
electrode position and size. This result was also confirmed in our numerical simulations. 

Once the optimal value of /3 is identified, synchronized voltage measurements and image capturing of the flapping 
plate are carried out using an oscilloscope and a high speed camera placed over the wind tunnel. Finally, the efficiency 
of the system is computed as Pe/Pf (see section [n|). 

Figure shows the evolution of the harvesting efficiency of the piezoelectric flag with the electrode’s length Lp 
for both experiments 1 and 2. A good qualitative agreement is observed between the experimental results and 
the numerical simulations corresponding to the same non-dimensional parameters. For experiment 1 (left panel), 
the maximum efficiency is achieved when the electrode occupies the whole length of the plate. For experiment 2 
(right panel), the maximum efficiency is achieved when the upstream extremity of the patches is located at s « 0.1. 
Quantitative discrepancies between experimental and simulation results are, as mentioned previously, due to errors in 
the estimation of the plate rigidity B used in our simulation. Although an effort has been made during the experiments 
in order to reduce the impact of the electrode sectioning on the plate rigidity, this effect exists, and it is not accounted 
for in the theoretical model of the plate. Discrepancies in the absolute efficiency values between the two experiments 
have their origin in the difference between coupling coefficients of each experiment: a = 0.073 for experiment 1 and 
a = 0.062 for experiment 2. 

The good agreement between experiments and simulations also validates the model presented in section and 
motivates its use in the parametric study of more complex geometries presented in the remainder of this section. 


B. Parametric study 


In this subsection, simulation results for the energy harvesting efficiency are presented in the three cases sketched in 
figure 

I: One single electrode with variable length Lp and variable position along the flag surface (two-parameter problem) 
II: Two electrodes of variable length covering the whole surface of the flag (one-parameter problem) 

III: Three electrodes of variable length covering the whole surface of the flag (two-parameter problem) 


For each case, we analyze the influence of size and distribution of the electrodes on the energy harvesting efficiency, 
defined as in the experimental case and expressed in non-dimensional form in Equation (24). 
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FIG. 6. Three cases of study in the energy harvesting efficiency of the piezoelectric flag. 





FIG. 7. Energy harvesting efficiency of a single electrode with Lp = 1 as a function of P for M* = 0.6, U* = 13 and a = 0.25 
(solid fine) and M* = 10, U* = 8.25 and a = 0.25 (dashed line). 


A strong dependence of the efficiency on the system’s physical parameters such as the mass ratio M*, /3, U* and a 
is observed m- Before exhibiting results of geometry optimizations, it is necessary to select a representative set of 
these parameters. First of all, due to the weakly-nonlinear nature of the model used in the present work, the accuracy 
of the results is restricted to flow velocities U* near the flapping instability threshold. Far from threshold this model 
exhibits discrepancies with previous simulations and a fully-nonlinear model of the system should be considered m- 
Thus, all the results presented further in this work correspond to near threshold situations {U* is close to U*), and 
hence, depend on the value of M* selected. 

Secondly, two typical values of the mass ratio M* will be considered: M* = 0.6, which corresponds to a plate in 
a low density fluid (i.e. flag in air), and M* = 10 which corresponds to a plate in a dense fluid {i.e. light plate in 
water). For each case, the flow velocity U* is adjusted to a typical value just above the instability threshold: U* = 13 
and 8.25, respectively. 

As explained in the previous section, maximum energy harvesting is achieved when the flapping timescale and the 
RC circuit timescale compare. Hence, the optimal value of /3 will depend on the flapping frequency, which in turn 
depends on the mass ratio M* [H]. On figure the efficiency 77 in a typical single electrode configuration is plotted as 
function of /3 for two different values of M*. One then observes that higher values of M* tend to reduce the optimal 
value of p. Indeed, while for M = 0.6 the optimal tuning parameter is /3 ^ 1 (as in the experimental case), for 
M* = 10 the value of /3 is about 0.3. In the following, these optimal values for P are retained. 

The influence of the coupling coefficient a, the electrode position and geometry on the flapping frequency has then 
been tested. No significant impact on the flapping frequency, and consequently on the optimal value of P has been 
observed. This is expected to result from the relatively low coupling coefficients a considered. 
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FIG. 8. Colormaps of the efficiency rj in the {Lp, 5) plane; for all cases, M* = 0.67, t/*=13 and /3 = 1. 


Case I : single electrode 

We consider here the case of one electrode of length Lp that partially covers the surface of the plate. It is attached 
at a distance S from the clamped end. Figure [^presents the efficiency of the harvesting system as a function of S and 
Lp, for hxed values of M* = 0.67, U* = 13 and /? = 1, and increasing values of the electro-mechanical coupling factor 
a. Each colored rectangle of these plots corresponds to one single numerical simulation. 

We observe that regardless of the value of a, the maximum efficiency is found for values of <5 7 ^ 0 and electrode 
lengths 0.35 < Lp < 0.5. In these cases, the electrode is mainly located on the downstream half of the clamped-free 
plate. Moreover, the maximum efficiency region is observed to change as the electro-mechanical coupling factor a 
is increased. For small values of a 0.1), the maximum efficiency region on the {Lp,S) space extends from long 
electrodes Lp = 0.8 placed at <5 = 0.2 to smaller electrodes of length Lp ~ 0.4 and S ^ 0.45 (Figure]^). As the 
coupling factor a is increased the maximum efficiency is also increased. Simultaneously, the maximum efficiency 
region shrinks and migrates in the {Lp,5) plane to a new optimal configuration for an electrode of length Lp ~ 0.4 


























12 



_ 0.08 
1 0.075 

-- 0.07 
-- 0.065 
-- 0.06 
-- 0.055 
-- 0.05 
- - 0.045 
-- 0.04 
-I 0.035 




FIG. 9. Normalized flapping amplitude difference AAn = {A — Aq)/Aq in the {Lp,5) plane for (a) a = 0.1 and (b) a — 0.5 
respectively, where Aq corresponds to the flapping amplitude with no coupling (a = 0). (c) Normalized critical fluid velocity 
AUiq = {Uc — Uco)/Uco as a function of the electrode position in the {Lp, S) plane where Uco corresponds to the flapping critical 
flow velocity with no coupling (a = 0). (d) A9i integrated over one oscillation cycle for a plate with no piezo-electrical coupling. 
The mode shape of the fluttering flag over one oscillation cycle is depicted in the inset. Other parameters are M* = 0.67, 
f/* = 13 and /? = 1. 


and S ^ 0.55 (Figure]^). 

To explain this difference in the optimal configuration, the influence of the geometrical parameters (d, Lp) on the 
flapping amplitude is represented on figure [^a,b) for small and large values of a. These figures show the normalized 
relative flapping amplitude AA^ = (A — A(f)/AQ in the {S,Lp) plane for a = 0.1 and a = 0.5, respectively, where 
Aq is the flapping amplitude obtained with no coupling (a = 0). Comparison of these figures show that a only 
impacts the general magnitude of the amplitude fluctuation, and does not modify the parameter values leading to 
maximum (or minimum) amplitude. In particular, the maximum flapping amplitude corresponds to the optimal 
efficiency region achieved for high coupling. At large coupling, maximizing the flapping amplitude seems therefore to 
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be the optimal strategy. These maximum flapping amplitudes are also associated with the largest reduction of the 
critical flow velocity (Figure]^) which is consistent with the generally observed result that the flag flapping amplitude 
is a monotonically increasing function of U* — U* in the vicinity of the threshold. 

On the other hand, when the coupling coefficient is small, the flapping flag dynamics is only marginally modified by 
the energy transfers to the piezoelectric patches. In that case, and neglecting the feedback coupling of the piezoelectric, 
the amplitude of the equivalent current generator (Figure [^) can be directly obtained using Equation (§ from the 

kinematics of the flag in the case of no coupling. More specifically it is proportional to the relative rotation rate 

of the two edges of the piezoelectric patch. Maximum efficiency is therefore achieved for low coupling when is 
maximum, where 

A^=(|0(s+)-0(s-)|), (25) 

a result consistent with Figures a) andj^d). 

Before presenting some results at a higher value of the mass ratio, let us summarize the results obtained so far: 

— For weakly coupled systems (a = 0.1), the piezoelectric coupling has no influence on the flapping amplitude and 
the efficiency scales with A9i (figure [ojl) . 

— For strongly coupled systems (a = 0.5), the piezoelectric coupling a has an influence on the flapping amplitude 
which in turn also influences the efficiency. It was finally shown that the effect of a on AAjv can be directly 
deduced from that of a on AUn. 

For M* = 10 and a = 0.25, we observe different local maxima of the efficiency in the {6,Lp) plane. They correspond 
to an electrode of length Lp ~ 0.2 at (5 ^ 0.3, 0.5 (absolute maximum) and 0.75 respectively. The presence of several 
local maxima with shorter electrodes is explained by the fact that a larger M* leads to a higher order for the mode 
selected by the instability [H] : each local maximum corresponds to the localization of the piezoelectric patches in the 
zones of maximum curvature. This is confirmed by the strong correlation of the maximum efficiency with A9i during 
the oscillating cycle (Figure [T^). 

It should also be pointed out that the maximum efficiency for M* = 10 is about 2.5 times bigger than that for 
M* = 0.67. This result stems from the difference in mode structures at higher mass ratio: a shorter wave length 
penalizes long electrodes that only react to the average flag curvature along their length. It also suggests that for higher 
M* a larger number of smaller electrodes occupying the remaining free space on the plate surface could eventually 
increase the total efficiency even more. This aspect is investigated further in the following. 


Case II: two electrodes 

The case of two electrodes covering the entire plate is now considered for the same values of M*, U* and /3 as in 
previous case, and depends on a single parameter namely the length of the first piezoelectric electrode Xc (see figure 
[^. For M* = 0.67, the maximum efficiency is obtained for Xc ~ 0.35 (Figure [TT| a) and is approximately twice as 
large as the efficiency obtained for a single electrode covering the entire flag {xc = 0 or ccc = 1), and represents a 
6% increase of the efficiency in comparison with the optimized one-electrode configuration detailed in Case I. In this 
new optimal situation, the second electrode covers the most significant part of the plate’s surface. Also, it is observed 
that for the same electrode’s length, the downstream electrode leads to the largest energy harvesting, a result that is 
somehow expected because this region corresponds to the largest deflections of the flag. 

For M* = 10, two local maxima are observed around at Xc ~ 0.5 and Xc ~ 0.7. Figure confirms that the 

downstream piezoelectric pair contributes the most to the efficiency for each maximum. However, the maximum 
efficiency of the two electrodes configuration is slightly smaller than that of case I, a consequence of the full coverage 
of the plate by the two electrodes, in particular near the leading edge where the local rotation vanishes. Due to the 
dominant mode structure at large M*, shorter electrodes (either when using three or more electrodes or by choosing 
to cover only a fraction of the plate’s surface) are therefore expected to improve again the efficiency. 


Case III: three electrodes 

We finally turn to the case of full coverage by three electrodes, described by the positions of the left edge of the first 
and second piezoelectric patches, Xci and Xc 2 respectively. 
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(c)Aei (a = 0) 

FIG. 10. (a) Evolution of the efficiency rj with {Lp,5) for a = 0.25, M* = 10, U* = 8.25 and H* = 0.5. (b) Normalized 

critical fluid velocity AU^ = {U^ — U*q)/Uc- (c) A6i integrated over one oscillation cycle for a plate with no piezo-electrical 
coupling. The mode shape of the fluttering flag over one oscillation cycle is depicted in the inset. 


For M* = 0.67 the maximum efficiency is reached for Xci ~ 0.4 and Xc 2 ~ 0.9 (left plot in figure 12). In this 
configuration, the second electrode is approximately the same as in the maximum efficiency configuration of case I. 
The presence of two additional electrodes results however in a 20% increase of the efficiency in comparison with case 

I. 


For M* = 10, the arrangement of maximum efficiency {xci = 0.5,a;c2 = 0.7) is again related to the optimal 
configuration of case I, with the second electrode of length Lp = 0.2 positioned at 5 = 0.5. Taking advantage of the 
deformation mode, the third electrode of length Lp = 0.3 and positioned at (5 = 0.7 contributes in almost the same 
amount to the efficiency as the second electrode. As a result, the total efficiency of the three electrode configuration 
is about three times higher than that of cases I and II. 

This clearly demonstrates that for higher mass ratios, a larger number of electrodes results in significant improve¬ 
ments of the efficiency, while this strategy only leads to moderate improvements for lower M*. 
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xlO-4 



xlO-4 



Xc 

(a)M* = 0.67 


Xc 

(b)M* = 10 


FIG. 11. Energy harvesting efficiency for a piezoelectric plate covered by two electrodes of variables size for (a) M* = 0.67 and 
U* = 13, and (b) M* — 10 and U* = 8.25. In both cases a = 0.25 and H* = 0.5. Each figure shows the total efficiency of the 
two electrode configuration (solid), as well as the efficiency of electrode 1 (dashed) and electrode 2 (dash-dotted). 



(a)r; (M* = 0.67) 


(b)r; (M* = 10) 


FIG. 12. Energy harvesting efficiency for a piezoelectric plate covered by three electrodes of variables size, (a) : M* = 0.67 
and U* = 13. (b) : M* = 10 and U* = 8.25. In both cases a = 0.25 and H = 0.5. 


V. CONCLUSIONS 


This article investigated the possibility of converting kinetic energy of a flow into electrical energy through elastic 
deformation of a fluttering plate using piezoelectric effect, a problem that has received growing attention in recent 
years. A critical question lies in the optimization of the energy transfers between the fluid, solid and electric circuit. 
The present work focused on the role and influence of the location of the piezoelectric electrodes on the efficiency of 
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the harvester, explicitly describing the discrete nature of the piezoelectric coverage of the flag and hence, allowing for 
a better representation of the experimental and practical situations. 

To address this problem, the present work used a combination of theoretical analysis and experiments. A weakly 
nonlinear model considering a plate with a discrete number of piezoelectric patches pairs was developed and confronted 
to the experimental results. It was then used to determine geometries that maximize the energy harvesting efficiency 
in three different test cases: (1) one electrode partially covering the plate, (2) two electrodes fully covering the plate 
and (3) three electrodes fully covering the plate. For each configuration, the optimal geometries were described and 
analyzed. 

The results presented in this work emphasize the critical role of the piezoelectric elements positioning: because 
the electric charge transfer responsible for the powering of the output circuit is solely given by the relative change 
of orientation of the flag between both ends of the patch, the placement of this patch must be carefully designed to 
take full advantage of the non-uniform flag deformation and its specificities, in particular the dependence of the mode 
shape with the mass ratio. In the case of a single electrode and small values of the mass ratio, the present study in fact 
shows that a greater amount of energy can be harvested when positioning the electrodes on the downstream half of the 
flag. At a higher value of the mass ratio, more than one optimal positioning are found, which is a consequence of more 
spatial oscillations of the unstable mode shape. Using a larger number of electrodes allows to improve the efficiency 
of the system particularly for higher mass ratios, where the modal structure corresponds to a shorter wavelength and 
is therefore better suited to shorter electrodes that are able to capture the flag’s deformation best. 

The work in this paper extends the results obtained in the context of damping or energy harvesting of vibrating 
plates in still fluid [33l 02] to flow-induced instabilities and more than one piezoelectric patch. The present results 
allow to conclude that the optimal location for piezoelectric energy harvesting along the plate is on the downstream 
half, and in that regard are reminiscent of the conclusions of Ref. [S] obtained for a pure damping model of the energy 
harvesting process. In both cases, the local deformation of the flag is driving the energy harvesting process, either 
through changes in the local curvature or in the relative orientation at both ends of the piezoelectric patches, and 
the response of the deformation of the flag to the presence of energy harvesting plays a critical role in determining 
the optimal location for large coupling coefficient. Both approaches however differ in the relevant measure of flag 
deformation (either curvature or relative orientation), which also impacts the detail of the optimal location for the 
damping or piezoelectric patches. 

In the case of weakly coupled systems, it was also shown that the amount of harvested energy can be directly 
predicted from the kinematics of the flapping flag obtained for zero piezoelectric coupling. Maximizing the harvested 
energy hence requires to determine the positioning that will ensure a maximum relative variation of the orientation 
between the two ends of the piezoelectric components. In the case of strongly coupled system, the positioning of the 
electrodes may influence the flapping amplitude, and consequently the efficiency, because of the feedback effect of 
the piezoelectric patches on the flag dynamics. Near the stability threshold, the flapping amplitude can be directly 
related to the critical velocity fluctuations: in that regard, linear stability analysis may be a great tool to estimate 
the impact of the piezoelectric coverage on the efficiency in the case of large coupling. 
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Appendix A 


The vertical displacement y is expanded on clamped-free beam eigenmodes (l>p{s), 


N 

2 = 1 


(Al) 


After truncation to N linear modes, equation (15) is projected along mode p. Projection of the forces L{y), /m, /b, 
f^, /res and /reac are given in this appendix. Projection of L{y) along mode p reads, 
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For fmiy) we have: 
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The term /b(2/) projects as: 
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The mechanical-electrical coupling term gives 
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For the resistive term free we have : 
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where 


ns N 

fc = / </’p X! + MiMjqj + 4>jqj)ds. 

.• .-I 
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Notice that unlike the other coefficients is a function of time. It is due to the presence of the absolute value. 
This coefficient must then be recalculated at each time step. 

Finally, for /reac we have: 
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